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1 Introduction 

The cross section for 77 -^ tt^tt'^ and for 77 -^ tt+tt" has been calculated some 
time ago fl], ^ in the framework of chiral perturbation theory (CHPT) |^-[0] and 
of dispersion relations. For charged pion-pair production, the chiral calculation |1| 
at next-to-leading order is in good agreement with the Mark II data ||] in the low- 
energy region. On the other hand, for 77 -^ tt^tt^, the one-loop prediction |]l], |^ 
disagrees with the Crystal Ball data p| and with dispersion theoretic calculations 



|T^-|16] even near threshold. 

In the process 77 -^ tt+tt", the leading contribution^ is generated by tree dia- 
grams. One has a control on higher order corrections in this case, in the sense that 
it is explicitly seen that the one-loop graphs do not modify the tree amplitude very 
strongly near threshold [|l|. Tree diagrams are absent for 77 — > vr^vr" which starts out 
with one-loop graphs. It is the aim of this article to establish the region of validity 
of the chiral representation of this process by evaluating the amplitude at two-loop 
order. 

Is a next-to-leading order calculation sufficient in this case? If the corrections are 
large, the reliability of the result is certainly doubtful. However, a glance at the data 
shows that the corrections needed to bring CHPT and experiment into agreement 
are not large-a 25-30% change in amplitude is sufficient. Corrections of this size are 
rather normal in reactions where pions in an isospin zero S-wave state are present 
||T9| . As an example we mention the isospin zero 5- wave vrvr scattering length, whose 



tree-level value |^ receives a 25% correction from one- loop graphs [Q. Corrections 
of a similar size are present in the scalar form factor of the pion [^ . 

The amplitude for 77 -^ tc^tc^ also describes Compton scattering on neutral pions 
by analyticity and crossing. Do sizeable corrections in 77 -^ tt^tt^ then also show 
up in 77r° —>■ 771° ? Since there are no strongly interacting particles in the final 
state in this case, one might be led to suspect that the one-loop amplitude is a good 
approximation for this reaction. We find it interesting that this is not the case-the 
corrections to the leading-order term are in fact very large in this channel. 

The electromagnetic polarizabilities characterize aspects of the inner structure 
of hadrons. With the two-loop expression for the amplitude at hand, it is straight- 



^ In this article, we denote the first nonvanishing contribution to any quantity by "the leading- 
order term" , independently of whether it starts out at tree level or at higher order in the chiral 
expansion. 



forward to evaluate the polarizabilities at next-to-leading order in the quark mass 
expansion. Renormalization group arguments show that this expansion contains log- 
arithmic singularities of the type M^- In M^ and M,r In Af,r, and an order of magnitude 
estimate reveals that these contributions may easily be as large as the leading-order 
term, unless the relevant Clebsch-Gordan coefficient is small. We find that the latter 
is the case. 



Recently, a reformulation of CHPT has been given |^, where the effective la- 
grangian includes into each order additional terms which in the standard CHPT 
(considered here) are relegated to higher orders. To all orders, the two perturbative 
schemes are identical-in each finite order, they may, however, substantially differ. 
For an analysis of the process 77 -^ ti^ti^ in this generalized framework we refer the 
reader to Ref. |]23 . 



The article is organized as follows. In section 2, we set up the notation. In section 
3 we describe the low-energy expansion in a general manner and outline the specific 
procedure for the two-loop case in sections 4 and 5. The low-energy constants which 
occur in the amplitude for 77 -^ vr^vr^ at two-loop order are determined in section 
6. Section 7 contains a discussion of the amplitude and of the cross section at two- 
loop order. The Compton amplitude and the pion polarizabilities are described in 
section 8, whereas section 9 is devoted to a comparison of the chiral expansion with 
the dispersive analysis of 77 —>■ 7r°7r° by Donoghue and Holstein [|13]. Finally, a 
summary and concluding remarks are presented in section 10. 

2 Kinematics 

The matrix element for pion production 

7(91)7(92) ^ vr°(pi)7r°(p2) (2.1) 

is given by 

< 7r°(pi)7r°(p2)out I 7(9i)7(92)in >= ^{2nyS\Pf - P,)T'' , (2.2) 

with 



T^ = e^e^e^V; 



fii^ ; 



V,, = I rfxe-*(^i^+''^^) < 7r°(pi)7r°(j92)out | Tj^(x)j,(i/) | > . (2.3) 



Here j^ is the electromagnetic current, and a = e'^/An ~ 1/137. The decomposition 
of the correlator V^u into Lorentz invariant amplitudes reads with q\ = q^ ~ ^ (^^^ 
appendix ^ 

s 

T2f,u = 2sA^Ai, - z/^^f^j, - 2u{qi^Af, - q2^,A^) , 

Ts^iy = qi^q2u , 

T4^,u = s{qi^Ay - q2vA^) - u{qi^qi^ + q2^^q2u) , 

A^ = (pi-P2)m , (2.4) 

where 

s = (gi + g2)^ t = (pi- gl)^ u = {p2- qif , 

z/ = t-u , (2.5) 

are the standard Mandelstam variables. The tensor V^y satisfies the Ward identities 

gfV = ?2 V = • (2.6) 

The amplitudes A and B are analytic functions of the variables s, t and m, symmetric 
under crossing (t,M) — ^ {u,t). The quantities C and D do not contribute to the 
process considered here (gauge invariance). 

It is useful to introduce in addition the helicity amplitudes 

H++ = A + 2{4M^-s)B , 

H,_ ^ ^SMl^^B. (2.7) 

s 

The helicity components -ff++ and Hj^^ correspond to photon helicity differences 
A = 0, 2, respectively. They have partial wave expansions involving even J > X, 



H++ - 


- >; 

J=0,2,4.. 


hiis)d^oo{cos 


H+- - 


= >; 


h'^{s)d2Q{cos 



, _., , , (2.8) 

J=2,4,6... 



where 6 is the scattering angle in the center-of-mass system, qi ■ pi = |^i||pi| cos^. 



With our normalization of states < pi | P2 >= '2{2TT)'^p'l6^{pi—p2), the differential 
cross section for unpolarized photons in the center-of-mass system is 

H{s,t) = I i7++ p + I i/+_ p , 

/?(.) = il-AM'Js)'/'. (2.9) 

The amplitude for Compton scattering 

7(gi)7r°(pi) -> 7(?2)7^°(P2) 

may be obtained by crossing. In the center-of-mass system, the cross section for 
unpolarized photons is 

with 

s = {qi+PiY , i= {q2-qif ■ 
Finally, the optical theorem in the Compton channel reads with our phase convention 

e^Im5|..o,., = ^^^fS{s) . (2.11) 



This relation fixes the phase of A through Eq. 

The physical region for the reactions 77 -^ -k^-k^ and 77r° -^ 771° is displayed 
in Fig. 1, where we also indicate with shaded lines the nearest singularities in the 
amplitudes A and B. These singularities are generated by two-pion intermediate 
states in the s, t and u channel. 

3 Low-energy expansion 

We consider QCD with two flavours in the isospin symmetry limit ?7i„ = tti^ = m 
and equip the underlying lagrangian with hermitean, colour neutral external fields 
v,a,s and p in the standard manner, 

^ = ^QCD + ^7''(^M + fl/xTs)? - q{s - il5P)q ■ (3.1) 

Here C^^.^^ denotes the QCD lagrangian at zero quark mass, whereas m is contained in 
the scalar field s{x). The lagrangian ( p.l|) is invariant under local SU{2)l x SU{2)r x 
U{1) transformations 

g^^[(l+75)^ij + (l-75)^L]g (3.2) 



with 

Qr^l = e'^VR,L , 
Vr,l e SU(2) , 

= diag(0o, 0o) , <^o e R , (3.3) 

provided that the external fields are subject to the gauge transformation 

r'f, = gRTf^gli + igRdf^gli , 

I'f, = gLif^gi + Widf^gi , 

s' + ip' = gRis + ip)gl , 

r^, = v^ + a^ , l^ = v^-a^ . (3.4) 

Since the charge is not a generator of SU{2), we consider in the following the case 

{a,} = , (v,) y^ , (3.5) 

where (A) denotes the trace of the matrix A. The condition ( p.5|) is consistent with 



the transformation law (|3.4|) . The Green functions of the theory are generated by 



the vacuum-to-vacuum amplitude 

iZ{v,a,s,p) _ n I n > (Qa) 

The generating functional Z admits an expansion in powers of the external momenta 
and of the quark masses 0-]^, 

^ = ^2 + ^4 + ^6 + ... , (3.7) 

where Zn denotes a term of order E"^. We write the corresponding expansion of the 
amplitudes as 

I = l2 + h + h + ... ; / = V', A, B , (3.8) 

where it is understood that Z^ generates /„|. To calculate V^'^ , we set 

s = rhl , Vf, = QVf, , p = T^p , a^ = , (3.9) 

where 

Q = ^diag(2,-1) (3.10) 

is the charge matrix, and where t;^ and p denote flavour neutral external flelds. V^'^ 
is obtained from the term of order v'^p'^ in Z. 

^Notice that In is not of order E". 



3.1 Terms at order E'^ 

In the meson sector, Z2 is given by the classical action 

Z2 = dxC2{U,v,a, s,p) , (3-11) 

where £2 is the nonlinear a-model lagrangian 

£2 = ^{D.UD^U^ + x^U + xU^), (3.12) 

evaluated at the solution to the classical equation of motion S J C2 = 0. The 2x2 
unitary matrix U contains the pion fields, 

U = a + t-,a+ — = l, 

It transforms as 

U ^ gnUgi (3.14) 

under G = SU{2)l x SU{2)r x f/(l). The covariant derivative is 

D^U = d^U-ir^U + iUl^ , (3.15) 

and the field x denotes the combination 

X = 2B{s + ip). (3.16) 

F is the pion decay constant in the chiral limit, F.,^ = F{1 + 0{rh)), Fj^ ~ 93 MeV, 
and B is related to the order parameter < | gg | >. The physical pion mass is 

M2 = M2(l + 0(m)) , 

M^ = 2mB. (3.17) 

£2 is referred to as the effective lagrangian at order E'^. 

The term of order 0{v'^p'^) in the classical action Z2 vanishes and, therefore, one 
has 

^2 = ^2 = Vf = . (3.18) 



3.2 Higher orders in the energy expansion 

At higher orders in the energy expansion, the effective lagrangian consists of a string 
of terms. Reintroducing momentarily ^, one has 

C,s = C2 + hCi + h'^Cii + ---. (3.19) 

Here £4 contains all possible contributions with four derivatives, or two derivatives 
and one field X; o^' X^i ^ind similarly for the higher order terms (£4 contains in 



addition the Wess-Zumino-Witten lagrangian £wzw O)- The generating functional 
is given by 

and its low-energy expansion is obtained from Z = Z2 + hZ^ + • • • . One expands 

Ci = Ci + Cii + ]^iDji + Eie + Fie + --- -1 = 2,4,... , (3.21) 

where £/ denotes the lagrangian Cj, evaluated at the solution to the classical equation 
of motion 6 J C2 = 0. (To simplify the notation, we have dropped the SU{2) - indices 
in ^ and in the operators Cj, Dj, . . . .) The fluctuation ^ is of order h^'"^. Then one 
obtains 

^ J [{E^e + nC,0x{E2e + fiC.Oy] dxdy 

+ ^ [ \2F2^^ + HD^^] dx + 0{h^) , (3.22) 

2n J '- J^ 



A = 1-— o 
2h 



with S^i = J dxC^ff. 

At order E'^, this result amounts to evaluating one-loop graphs generated by £2 
and adding the tree graphs from £2 + ^£4 @]- These contributions then add up to 
Z4, which contains the leading-order term Kf'^. It is a specific feature of the process 
77 -^ tt'^tt^ that the counterterms contained in £4 do not contribute to V^^^^-the sum 
of the one- loop graphs is therefore ultraviolet finite ||l], H . 

The diagrams which generate Zq are displayed in Fig. 2. The solid-dashed lines 
stand for the propagator D2^, and the framed symbols I denote vertices from Cj 
according to Eq. (|3.21|) . 

In order not to interrupt the argument, we relegate the discussion of the lead- 
ing contribution V^''^ to appendix ^ and continue in the following section with the 
evaluation of the next-to-leading order term Vq'^. 
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4 Renormalization 

The evaluation of Vq'^ is complex. We outline in this and in the following two sections 
the procedure - omitting, however, all details. 

4.1 The lagrangians £4 and Cq 

The lagrangian £4 contributes to Vq"^ through one- loop diagrams, see Fig. 2. Its 
general form is M 



£(^) = J2kPi + --- , (4.23) 



where 



i=l 



Pi = \{u'u,r , 

P2 = ^{u^u,){u''u'') , 

Pi = ■^{Uf.X-) , 

Pg = ■^{f+[u^,u^]) , 

Pr = -^{X~r. (4.24) 



Here we used the notation 



u^ = iv) DfjUu^ = —iuD^U'^u = v)^ , 



X± = u'^xu^^ux'u , 

/r = uFru^±u^F^''u , (4.25) 



with u^ = U. The quantity Fj^'^ {F'^") stands for the field strength associated with 



the nonabelian external field v^ + a^ (v^ — a^) 



The ellipsis in ( [4.23| ) denotes polynomials in the external fields which are inde- 



pendent of the pion variables. These do not contribute to S'-matrix elements and are 



therefore not needed in the following. Finally, the anomaly term £wzw contributes to 
Vg'^ mi- This is beyond the accuracy of the low-energy expansion considered here. 



The realization of G on m is 



u 



.V /-( — n I h\ I ht n I /~t ' 



gRuh^ = hugi , (4.26) 



such that 

/ ^ hlh^ (4.27) 

for the quantities in ( [4.25| ). The low-energy constants /j are divergent, except l^. They 
remove the ultraviolet divergences generated by the one-loop graphs-we discuss them 
in more detail below. 

In the construction oi C^^\ the equation of motion 5 J C2 = has been used. It can 
be shown that adding terms to C^'^^ which vanish upon use of the equation of motion 
affects the generating functional at order E^ by a local termj^these contributions 
may thus be omitted. 

The lagrangian Cq contributes a polynomial part to V^'^ which cancels the ultra- 
violet singularities generated by the two-loop diagrams. The general structure of Cq 
is not yet available in the literature [^ . Concerning the present calculation, we note 
that the lagrangian 

^e = T^\f+iJ.pf+ + f-tipf- )T„ + ■ ■ ■ , 
Tp^ = di{upU„) + gp„{d2{u^Up) + d'i{x+)} (4-28) 

generates a polynomial in Ag, Bq which has the same structure as the divergent part 
in the two-loop contribution, 

20 

^6 = ^[16(^3 -t/2)M2 + (rfi + 84)s] + --- , 

B, = -^A + ---- (4-29) 

We may therefore remove the divergences in V^''^ by simply dropping the singular 
parts in Aq and Bq, see below. 



4.2 Regularization and renormalization 

We use dimensional regularization and set 

uj = d-4: , (4.30) 



•^We thank G. Ecker for an explicit proof of this statement and for ilhiniinating discussions 
concerning the material in section 4.2. 
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where d is the dimension of space-time. We introduce the renormahzation scale 
/i such that the calculation is scale-independent at each step. In the following we 
outline the procedure P7| . 

Consider the couplings U in iZ'^^) which carry dimension (mass)"^. We treat U as 
/^-independent parameters by writing in the minimal subtraction scheme 

/, = /.- 1^ + r + ^iTi + 0{uj')\ ; z = 1, . . . , 7 , (4.31) 

with 

J/MS JIMS 

/^^ = -5^ , ^^-^ = -r • (4.32) 

The divergent terms 6i/uj remove the one-loop singularities in Z4 (the 6i are related 
to the 7j used in g] by (5j = •yi/lGrr^). Since the k occur in Zq via loop insertions, 
the constants /^^ in general also contribute at two-loop order. For an illustration of 
the renormahzation procedure, we consider the amplitude F'^B which has dimension 
(mass)^'^. We write for the contribution from the loops 

B = ^ {Laurent-series of (fi'^'^F^B) at cj = o} . (4.33) 

From C'^'^\ only I2 contributes, 

ij {M'-f{s/M^ t/M^- u) + M-/2 g{s/M\ t/M'; cu)} , (4.34) 

where / and g are singular as cu — > 0, 

/ = ^ + - + fo + Oiiu) , 

C(J^ UJ 

9 = - + 9o + ujg+i + 0iiu^). (4.35) 

00 

The Laurent-series ([4.33|) becomes 



« p. 



g^^J4 + ^''" + '^" +,8"' + OM|. (4,36) 



The residues of the pole terms are 



/?2 = ^ , A = ^1 , A,i = '^2^0 + /i , (4.37) 



and for the finite part we find 



^Ms^^ J^^l^ ^4;ms 1^ ^^^1^ ^^^ . (4.38) 

8 [ /x^ J /^ A* 
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Here the q stand for linear combinations of /i, . . . , g+i. 

Without having evaluated any Feynman diagram, we have already obtained sig- 
nificant information on the structure of the two- loop result Q: 



1. The residues of the pole-terms in (|4.36| ) are polynomials in the external mo- 



menta and in the masses on general grounds [^. For dimensional reasons, 
these polynomials reduce to pure numbers in the present case. In addition, 
from the cancellation of the logarithmic terms which are generated by expand- 



ing the factors M^^ and M'^ in ( [4.34 ), it follows that the residue of the double 



pole in / is proportional to the residue of the single pole in g, 

2/2 + 5291 = , (4.39) 

or 

2/?2 = S2(3i . (4.40) 

2. The amplitude / contains a nonlocal singularity fi/uj which is generated by 
divergent subgraphs. This nonlocality must be cancelled by the nonlocal sin- 
gular part 52g{)/oJ in the graphs generated by loops with C'^'^\ in such a way 
that /3i 1 = 8290 + /i becomes a pure number. 



3. As is seen from Eq. (|4.38|) , the finite part [3^^ contains chiral logarithms In M^ 
and InM^. The coefficient of the leading term In M^ is proportional to the 
residue of the single pole in g, whereas the linear piece InM^ is multiplied with 
a nonlocal function. These singular terms cancel the chiral logarithms in C2 and 
thus generate a smooth behaviour of the amplitude in the chiral limit m — i> 
(at fixed s, t 7^ 0). 

4. At the Compton threshold s = 0, t = M^, the quantities ci and C2 are indepen- 
dent of the quark mass. Therefore, the chiral logarithms In^ M^ and In M^ in 
Eq. ( [4.38|) remain, and we conclude that the finite part /5^^|s=o,i=M2 blows up 
in the chiral limit. In other words, the slope of the form factor Vg = 2sB (see 
appendix A) is infrared singular. 



^^U.^M. = ^ -52ln^+4/2- ln- + c,ln- 



where the Cj are independent of the quark mass. Notice that In^ M^ occurs 
together with l^^ In M^ in a particular combination which is dictated by Eq. 
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We now define the renormalized amplitude B""'^ as 

5- = i_(^Ms ^ ^MS) , (4.42) 

where the scale-dependence of b^^ is chosen such that B""'^ is independent of /i, 

^A^-n = . (4.43) 

(The low-energy parameter b^^ is the sum of the finite pieces of the relevant coun- 
terterms at order E^ in the effective action.) 

We have formulated the renormalization procedure in the minimal subtraction 
scheme, where ln(47r) and F (1) occur. One may eliminate these terms in the standard 
manner [|2U[. Below we use the conventions of [Q. The final result for Bq contains 



one unknown new parameter b^, 

Analogously, the renormalized amplitude Aq contains two unknown new parameters 
al and a'2. 



The ellipses in ( [4.44| ) and ( ^.45| ) stand for the finite contributions from the loop- 
integrals. 

5 Evaluation of diagrams 

Here we discuss further aspects of the two-loop calculation. 

5.1 The diagrams 

It is straightforward to generate from figure 2 the diagrams for 77 — > vr^Tr" at two- 
loop order-one has simply to insert photon and pion vertices in all possible ways. 
For illustration, we display one class of graphs in Fig. 3. The solid (dashed) line 
denotes charged (neutral) pions. The four-point function on the right-hand side is 
the d-dimensional elastic vrvr scattering amplitude at one-loop accuracy, with two 
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pions off-shell (the one-loop graphs for 77 —>■ tt^tt^ are thus also included in Fig. 3). 
The symbol d^l stands for integration over internal momenta with weight 

(5.1) 



[M2-(/ + giP][M2-(g2-/)2] ' 
where M^ denotes the physical pion mass in one-loop approximation, 



Mf = M^ 



A/f2 1 M2 



(5,2) 



The momenta of the charged pions running in the loop are (/ + gi)^ and (g2 — O'^- 
We do not display the remaining diagrams. 

5.2 Numerical evaluation of diagrams 

The derivative nature of the interaction makes the algebraic part of the calculation 
tedious. As for the numerical part, we have to evaluate the amplitudes in the phys- 
ical region for 77 -^ tt'^tt^ where branch-points and cuts appear. These render the 
numerical evaluation of the Feynman integrals nontrivial. To illustrate our proce- 
dure to cope with this difficulty, we consider the two-loop box diagram Fig. 4a. We 
write for the subdiagram (enclosed by dash-dotted lines) a d- dimensional spectral 
representation 

^ ^ i4M2 a-t ^ ^ 

where i contains loop-momenta. This leads to a spectral representation for the full 
diagram, 

PCX) 

dap{(7,d)b{s,t,a,d), (5.4) 



f 



where 6(s, t, a, d) denotes a one-loop (i-dimensional box diagram (one of the internal 
lines carries mass ^/ff). After removing the subdivergence generated by J{i,d), we 
obtain the finite part by writing a fixed-t dispersion relation. This procedure allows 
one to evaluate numerically the amplitude also in the region s > AM^,t < 9M^. 

5.3 Checks on the calculation 

As is shown in appendix ^, the most general expression for the amplitude e^eg V^^ 
contains five form factors which are linearly related through two Ward identities and 
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through Bose symmetry, 

Li{Vi) = 2Vo + sV4 - {t - u)V5 = , 

L^iV,) = sVy-{t-u)Vg = , (5.5) 

LsiVi) = V, + Vr = . 

The amphtudes A and B may be obtained from V4 and Vg, 

A = -14 , B = Vg/2s . (5.6) 

To get an optimal control of the calculation, we have fully evaluated all five form 
factors in d dimensions and have made the following consistency checks of the results: 

1. We have verified that the relations (|5.5|) are satisfied numerically within ma- 
chine accuracy in low dimensions (below threshold). 

2. We have determined the divergence structure at a; — »• 0, 



p{2) p(l) 



V, = ij'-{^ + ^ + R, + Oiuj)} , (5.7) 

and have verified that the residues Pj are polynomials in the external mo- 
menta and in the masses. These polynomials obey Eq. (|5.5| ) as well, 

L„(pf))=0 ; m = 1,2,3 ; fc = 1, 2 , (5.8) 



and are related in the manner discussed in the previous section, see Eq. ([4.40|) 



for P^^^ and P^'^\ 
3. We have checked that also the finite parts Ri satisfy numerically 

Lm{R^)=0 ; m = 1,2,3 (5.9) 

within machine accuracy below threshold. 

At this stage, we have written fixed-t dispersion relations for some of the finite 
parts in the manner mentioned above. This allows one to evaluate the complete 
amplitude in the physical region for pion-pair production. 



4. We have then worked out the S'-wave projection h^{s) ( p.8|) and have verified 
numerically that this amplitude has the correct phase at s > 4M^, given by the 
(tree plus one-loop) elastic vrvr scattering S'-wave phase shift (in the appropriate 
isospin decomposition). 

In addition, we have made many other cross-checks. 
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6 Low-energy constants at order E^ 

Once the program described above is carried through, one ends up with ultravio- 
let finite and scale-independent amplitudes A and B which contain the parameters 
F,Mt^; 11, III (* — 1)2,3,5,6); a[,a2 and 6^. F is related to the physical pion decay 
constant F^ [ffl, 



F^ = F 



A/f2 1 M^ 



(6.10) 



We may therefore replace F by F^^ at the expense of introducing l\. The expressions 
for the loop- amplitudes simplify if one uses @] the scale-independent parameters /«, 



2 ^'^ + ^"7 



/[ = ^(!, + ln— ). (6.11) 



Their values are displayed in column 2 of table 1. We note that l^ and Iq in the 
present application always appear in the combination 

Ia = h-h = 2.7 . (6.12) 

[The constant /a is related to the low-energy couplings Lg and L\q which occur in 
the SU{?>)l X SU{?>)r version of the one-loop amplitude 77 -^ tt+tt^ II3] t>y ^a = 
1927r^(Lg + L\q)] Next we observe that we may absorb /[j^'s into the low-energy 
constants at order E^, because they contribute a polynomial piece only. We are 
therefore left with a\, a\ and V as the only new unknowns. We estimate these in the 
standard manner 0, ^ , replacing them at a scale \i = 500 MeV • • • 1 GeV by the 
contribution from resonance exchange. Let 

r (/x) = E ^"^ + ^'(z") ; I = ai,a2,b , (6.13) 

where the sum denotes contributions from scalar, (axial-)vector and tensor exchange. 
Our estimate for F{Mp) consists in setting F{Mp) = 0. 

The quantities I^ are evaluated in appendix 0. The results of the calculation 
are displayed in table 2, where the individual resonance contributions I^ are listed. 
Column 6 contains the sums of those contributions which have a definite sign. 

To estimate the effects of the systematic uncertainties in the values of these 
couplings, it is useful to furthermore consider the helicity amplitudes -ff++ and -ff+- 
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Table 1: Phenomenological values [^ |31| and source for the renormalized coupling 
constants li ,i = 1 . . .6. The quantities 6i in the fourth column determine the scale 
dependence of the /[(/i) according to Eq. ( |6.11[ ). In the text we also use Ia = h~h, 
see Eq. ( |02|) . 



i 


k 


source 


167r2(5i 


1 


-0.8 ±1.2 


Ke4, VTTT -^ TTTT 


1/3 


2 


5.8 ±0.7 


Ke4,, TTTT -^ 7171 


2/3 


3 


2.9 ±2.4 


SU(3) mass formulae 


-1/2 


4 


4.3 ±0.9 


Fk/F^ 


2 


5 


13.8 ±1.3 


71 — > eZ/7 


-1/6 


6 


16.5 ±1.1 


<P >^ 


-1/3 



Table 2: Resonance contributions to the coupling constants a^, aj and If . Column 6 
contains the sums of those contributions which have a definite sign. The calculation 
is presented in appendix 0. 





jR 




jR 


r 


UJ 


P° 





A(l+-) 


Y.rI'^ 


S(0++) 


/2 


\f 


-33.2 

12.5 
2.1 


-6.1 
2.3 
0.4 


-0.1 
~0 
~0 


0.0 

-1.3 

0.7 


-39 

13 

3 


±0.8 

±1.3 

0.0 


±4.1 
±1.0 
±0.5 



and the corresponding low-energy constants /i!j_ and h 



rr2ioops 



rr21oops 



(T6^K^^ + ^^4 + 



\{M^-tu)_^^ ^ 



s(167r2F2)2 

a\ + 86" , /i^ = a^ - 26" , K 



(6.14) 



From column 6 in table 2 we obtain the central values of these couplings. According 
to experience with resonance saturation at order E'^, we associate a 30% uncertainty 
to the contributions generated by (axial-) vector exchange and a 100% error to the 
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contributions from scalars and from /2. Adding these errors in quadrature, we find 

= -14±5 , 



h\{M,) 



7±3 , 
3±1 . 



(6.15) 



Notice that tensor exchange does contribute neither to /i^ nor to /i^, because the 
couphng ( |D.^ ) is purely D-wave. Scalars do not affect h^_. 



In Ref. |^5|, these couplings have been determined i) from vector- meson exchange 
and using nonet-symmetry, and ii) from the chiral quark model, with the result 



{hi,h:,hL)\,=M, 



(—18, 9, 2) vector-mesons (nonet) 
(—12, 6, 2) chiral quark model 



(6.16) 



which agrees within the uncertainties with the values in ( |6.15| ). 



This completes the determination of the parameters which occur at two-loop order 



in 77 — > 7r°7r°. 



7 Amplitudes and cross section to two loops 

7.1 The amplitudes: analytic results 

We obtain the following expression for the amplitude A to two loops, 

A = A4 + Ae + 0{E^) , (7.1) 



or 



A = ^^{s-M'^) + Ua + P^ + 0{E'). 



(7.2) 



The unitary part Ua contains s, t and u-channel cuts, and Pa is a linear polynomial 
in s. Explicitly, 



U, 



2 ^ 



Gis) {s'-Ml)J{s) + C{s,k) 



ik - 5/6) 



+ 



// 



[s-M'^)J{s) 



+AA{s,t,u) , 



247r2F4 
sG{s) + 2Ml{G (s) - 3 J (s)) 



"00 



(7.3) 
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with 



C{s,h) = -^{2{h-4/3){s-2M'^f + {k-5/6){4s'-8sM^^ + 16Mt)/3 



-3M% + 12MXs - Ml)h - l2sMl + 15M^} 



"oo 



i(3cos0^-i; 



(7.4) 



The loop-functions J etc. are displayed in appendix |C[ The quantity G^(s) in 
Eq. (|7.2|) stands for G(s), evaluated with the physical pion mass, and 6 denotes the 
scattering angle in the center-of-mass system. The term proportional to rfgo i^ ^a 
contributes to D-waves only. For A^ see below. 

The polynomial part is 
1 



Pa 

ai 
a2 



{IGn^F^y 



;[aiM^ + a2s] , 



al + ^ {4/2 + l{8l2 + 12[a - ^) - Y^2 + 12h + ^^° 



9 



1 r - - 9 ^- - fiQ7^ 

— \P + l(2l2 + 121^ + -)- % + 12rA + ^^ \ 
18 I ^ 3^3 144] 



In 






The result for B is 



or 



with unitary part 



U 



B 



B = Ub + Pb + 0{E^ 
{h - 5/Q)H{s) 



2887r2Fis 



+ AB{s,t,u) . 



(7.5) 

(7.6) 

(7.7) 

(7.8) 



For the polynomial we obtain 
Pb 
b 



;i67r2F^2)2 



6" 

36 



h 



i'+i{2h + -)-i + 



393 



144 



(7.9) 



The integrals AA,B{s,t,u) contain contributions from the two-loop box and acnode 
diagrams displayed in Fig. 4. It turns out that these contributions are very small 
for the cross sections below ^/s < 400 MeV, both for 77 -^ 7r°7r° (0.1% at 400 MeV) 
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and for the crossed channel 77?" -^ 'yir^ (1.5% at 400 MeV). Therefore, one obtains 
a rather compact and accurate representation of the two-loop amplitudes by simply 
setting Aa,b = in Ua,b^- 

7.2 The cross section 77 -^ tt^tt^ 

In Fig. 5 we display the data for the cross section cr(s; |cos^| < Z = 0.8) as de- 
termined in the Crystal Ball experiment 0. They are shown as a function of the 
center-of-mass energy E = ^/s. The solid line denotes the two-loop result, evaluated^ 
with the amplitudes ( [7.2| )-( |T!9| ). For the low-energy constants a^ja^ and 6*" we have 
used the values from column 6 in table 2, and the values of k are the ones displayed in 
table 1. Shown is furthermore, with a dashed line, the one-loop result [Q], obtained 
by setting Ua,b = Pa,b = 0, see also appendix |B|. Finally, the dash-dotted hues 
display the result of a dispersive analysis (Fig. 23 in Ref. [0). In that calculation, 
use was made of the J = 0, 2 S'-wave tttt phase shifts from Ref. |H^ (these phase 
shifts satisfy constraints imposed by Roy- type equations [Q). 

The two-loop result thus agrees well with the data and with the dispersive analysis 
of Pennington |jl2| in the low-energy region. 

We found it interesting to see which contributions are responsible for the increase 
in amplitude and cross section near the threshold. In Fig. 6, we display with a solid 
line the cross section, evaluated at A^^^ = and without resonance exchange. [The 
change compared to the full result (solid line in Fig. 5) is 0.2 nh at E = 400 MeV 
and thus negligible.] The dashed hne corresponds to k = 0, and the dash-dotted line 
is obtained by setting h = 13 = 0. We conclude that the increase is due to I2, h and 
^A = h~h- To make this statement more quantitative, we note that the dependence 
of the cross section on the /j-values can been summarized with the expression 

a"°°p=(s) ~ A^or"°°p(s) , 

N = 1 + (-5.8 + 5.0ri + 4:.9l2 - 0.2/3 + 5.4r4 + 3.71a) ■ lO'^ 

= 1 - 0.058 - 0.040 + 0.283 - 0.005 + 0.232 + 0.100 

~ 1.51 (7.10) 

which is accurate to a few percent up to 450 MeV. The analogous expression for the 



"^An analogous result holds for the elastic tttt scattering amplitude, which contains unitarity 
contributions with t- and u-channcl cuts which are negligible below E ~ 500 MeV for S-waves. 
^We use Fyr = 93.2 MeV, M^ = M^o = 135 MeV, unless stated otherwise. 
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helicity amplitude i/++ at the physical threshold s = AM^ reads 

rr2loops AT" rriloop 

A^ = 1 + (2.5 + O-eli + 1.212 - 0.2r3 + 2.7/4 + 2.71^) ■ lO"^ 
= 1 + 0.025 - 0.005 + 0.068 - 0.006 + 0.114 + 0.073 
~ 1.27 , (7.11) 

with 

MlHlX" = 5.8 • 10~^ (7.12) 

The contributions from li,l2 are vrvr rescattering effects. They amount to a 24% 
increase in the cross section (out of 51%) and to 6% in i/++ (out of 27%). The 
renormalization of F^^ (contribution from I4) amounts to a 23% increase in a. 

7.3 The amplitudes: numerical results 

To get more insight into the characteristics of the two-loop corrections, we display in 
Fig. 7 the real and imaginary part of the helicity amplitudes ±10'^M^H^± at t = u. 
The solid line shows 10^M^if_|_+. It incorporates all contributions except Aa,b- The 
dashed line is the same amplitude for the one-loop case, and the dash-dotted line 
is for — 10^M^i7_|__ with the same input as for the solid line. The curves start at 
E = 2Mt^, and the crosses refer to the center-of-mass energy of the tt'^tt^ system in 
100 MeV steps. The amphtude H+^ changes very rapidly just above threshold and 
is nearly purely imaginary in the region 350 MeV < E < 400 MeV. As expected, the 
amplitude -ff+_ which starts out with a D-wave term is very small at low energies. 
Resonance exchange adds to i7++ a positive real part, thus increasing the cross 
section below ~ 400 MeV and decreasing it above this energy. 

In Fig. 8, we display the quantity 10^M^if++ at t = m as a function of s/M^. 
Above the threshold s = 4:M^, the modulus is shown. The solid (dashed) line denotes 
the full two- loop (one- loop) result. While the two-loop contribution to the modulus 
is below 30% in the threshold region, it modifies iJ++ substantially (percentage-wise) 
at s = 0, where the amplitude is small, see also the discussion below. Furthermore, 
we note that the zero at s = M^, which occurs in the one-loop approximation Eq. 
( p.6|) , is only slightly modified by the loop corrections. Finally, we display with o 
(+) the modulus of the S-wave projected part of -ff++, taken from Fig. 19 (23) in 
Ref. m. 
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7.4 Error estimates and range of validity of the chiral rep- 
resentation 

The uncertainty in the amphtude has two sources. Firstly, the low-energy constants 
li, /i!j_ and /i^ used above contain certain errors. For the /j, these are displayed in 
table 1. The systematic errors in the low-energy couplings at order E^ have been 
estimated in the previous section, see Eq. (|6.15| ). Secondly, we are concerned here 
with an expansion in powers of the quark masses and of the external momenta. 
Higher order terms in this expansion (three loops and beyond) will change the cross 
section accordingly. 

We discuss first the effect of the uncertainty in the low-energy constants and 
concentrate for simplicity on /;,!J_ and h^. In Fig. 9 we show the variation of the cross 
section according to the error estimates in Eq. ( |6.15|) . The calculation is done at 



Aa,b = 0. The dashed lines embrace the region generated by assigning all possible 
combinations of signs to the systematic errors in the couplings /i!J_ and h^ according 
to Eq. (|6.15| ). The dash-dotted line corresponds to the central value in ( |6.15[) . 



It is clearly seen that, below 400 MeV, the uncertainties in /i!j_ and in h^ do 
not matter. (Since we estimate the couplings with resonance saturation, this is 
a reformulation of earlier findings p4|-|^.) Varying the scale at which resonance 
saturation is assumed between 500 MeV and 1 GeV results also in a negligible change 
in the cross section below E = 400 MeV. Beyond this energy, the uncertainty becomes 
more pronounced. Because the contribution from H^_ is tiny (see Fig. 7), only h'^_^ 
and hi really count. One might thus be tempted to extract these couplings from 
more accurate data in the range E = (400 — 600) MeV. This would be interesting, 
because /i!^ determines the difference of the electric and magnetic polarizabilities of 
the neutral pion at two-loop order, see below. However, in this energy range where 
s ~ (9 — 20)M^, hi is much more important than /i!j_. It will, therefore, be rather 
difficult to extract /i!j_ reliably in this manner JlSl. On the other hand, it would be 



interesting to perform a combined analysis of the two related processes 77 —>■ tt^tt^ 
and rj —>■ 7r°77 [^| in the framework of SU{3)l x SU{3)ji in order to obtain maximal 



information on the low-energy coupling constants which enter these amplitudes^. 

Turning now to the corrections from yet higher orders, we use the fact that 
^2ioops^^iioop ^ (i + e)2 with e ~ 0.25 and estimate a/a"""" ^ (l-e)-^. This amounts 
to a 15 — 20% uncertainty in the two-loop result below 400 MeV. At higher energies, 
^We thank J. Bijnens, M. Knecht and J. Stern for pointing this out to us. 



22 



the error in the cross section is more difficuh to assess. It may well turn out, however, 
that a more precise determination of the low-energy couplings leads to the conclusion 
that the chiral representation of the amplitude at the two-loop level is even valid up 
to E = (600 - 700) MeV in this channel. 

8 Compton scattering and pion polarizabilities 

The amplitudes A and B are analytic functions of s and t. At s < 0, they describe 
Compton scattering, 

7(gi)7r°(pi)^7(g2)vr°(p2). (8.1) 

We discuss this reaction in the present section, where we also work out the (neutral) 
pion polarizabilities at next-to-leading order. 

8.1 Compton scattering 

The cross section 77 — >■ vr^vr'^ receives a substantial correction near threshold due 
to 7TTT final-state interactions - which are absent in Compton scattering. Are then 
the two-loop contributions small in this channel? Fig. 8 shows that this is not the 
case: in the one-loop approximation, the amplitude -ff++ is one order of magnitude 
larger in the 77 -^ tt^tc^ channel than at Compton threshold. Therefore, even tiny 
corrections in 77 — ^ tt^tt^ may appear large in Compton scattering [^. In Fig. 
10 we display the cross section a"''^ ^'^^ as a function of the center-of-mass energy 
E^T,. The solid line shows the result of the two- loop calculation and the dashed line 
displays the one-loop approximation. They differ by one order of magnitude already 
near threshold. This is mainly due to the effect of the low-energy constant /i1 in 
i/+_ (omega-exchange in the language of resonance saturation |^). Putting if+- to 
zero results in a'^'^ ~*'^'^ = 0.7 nb at E'^tt = 350 MeV (dotted line). Purely two-loop 
effects however also change the cross section by roughly a factor of two at E^^, = 350 
MeV (dash-dotted line, evaluated at /i!j_ = h^ = 0). 

In summary, the Compton amphtude is tiny at leading order, and it is therefore 
rather unstable against the corrections generated by higher order terms. 
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8.2 Pion polarizabilities 

For a composite system it is customary to include the electric and magnetic polariz- 
abilities among the fundamental parameters-such as the electric charge, the magnetic 
moment and the mass-characterizing the low-energy limit of the coupling with the 
photon in the Compton amplitude. Hadrons are no exception, hence the theoretical 
description of their dynamics can be tested through the experimental determina- 
tion of the hadron polarizabilities |^. To set notation, we first consider Compton 



scattering for charged pions, 

liqi)^^iPi) - 7(g2)vr+(p2) , (8.2) 

in the laboratory system p^ = M^^. (In order to simplify the notation, we use the 
symbol M^ to denote both the charged and the neutral pion mass.) The electric (a-,,) 
and magnetic (P-,,) polarizabilities are obtained by expanding the Compton amplitude 
at threshold, 

rpC 



ei ■ 6*2* ( — - a^uJiUJ2 ) - Pn {qi X el) ■ (g2 x €2") + 



3) 



a^±P^ = -—H^^{s = 0,t = M^^) , (8.4) 



with gf = (ui, Qi). In terms of the helicity components ( |B.4D , one has 

where the bar denotes the amplitude with the Born term removed^. For neutral 
pions, one uses the analogous definition, 

a^o±p^o = ^H+^{0,M^) , (8.5) 



or, in terms of A and B, 



a , . _ . .2 



- Q. 

PttO = -^^-4|^=o,t=Af2- (8-6) 

Below we also use the notation 

{a±/3f = a^±p^ , 

{a±p)^ = a^o±f3^o . (8.7) 

An unsubtracted forward dispersion relation for the amplitude B gives with ( p.ll|) 

and analogously for the charged channel. 

"^ In H^_ first set t — Mj, then s ^ 0. We use the Condon-Shortley phase convention. 
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8.3 Data on pion polarizabilities 

There exist up to now two determinations of charged pion polarizabihties via mea- 
surement of the Compton amphtude. At Serpukhov [0, radiative pion-nucleus 
scattering n^A^ tt^'jA has been used. Here the incident pion scatters from a 
virtual photon in the Coulomb field of the nucleus. In the pion production process 



7P ~~^ 77r"'"n examined at the Lebedev Institute [42|, the incoming photon scatters 



from a virtual pion. Analyzing the data with the constraint [a + (3)'-' = give^ 



c J 13.6 ±2.8 [41] 



(a-(3)'^ = { LJ (8.9) 

^ ^ ^ 40 ± 24 ra . ^ ^ 



The Serpukhov data have been analyzed also relaxing the constraint (a + /5)*" = 0, 
with the result 

{a + pf = 1.4±3.1(stat.)±2.5(sys.) [43] , 

{a-pf = 15.6±6.4(stat.)±4.4(sys.) H . (8.10) 

Here we have converted the value quoted for /3^ into a number for {a — (3)^ , adding 
the errors in quadrature. 

Furthermore, also the process 77 -^ mr may be used to obtain information on the 
polarizabilities. Since in this case the amplitude at low energies is mainly sensitive 
to S-wave scattering, only {a — (3)'"'^ can be determined from the presently available 
1^, ^ data. In Ref. []TB|, unitarized S-wave amplitudes have beeen constructed. 



which contain {a — (3)^'^ as adjustable parameters. A simultaneous fit to Mark II 
and Crystal Ball data gives 

{a-(3f = 4.8 ±1.0 [16] , 

{a -(3)^ = -1.1 ±1.7 [|1|] , (8.11) 



where we have taken into account that the definition of the polarizabilities in [|I^ is 
Air larger than the one used here, see El], Eq. 1. 



The value ( ^ . 1 1| ) for [a — (3)'^ is consistent with Refs. |4^ , ^ within 1^ standard 
deviations, but not consistent with [^. The large relative error in {a — (3)^ reflects 
the fact that the threshold amplitude 77 -^ tt^tt*^ is quite insensitive to large relative 
changes at the Compton threshold, as we discussed above (see also Fig. 10 in [0). 
The determination of (a— Z?)*"'^ from 77 -^ nn furthermore suffers from uncertainties 

^We express the polarizabilities in units of lO^'^^fm' throughout. 
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which we find difficuh to estimate in the approach used by Kaloshin and Serebryakov 
||T6| , which does not provide a systematic way to control the inherent uncertainties in 
the model amplitude used to fit the data. It might be interesting to merge dispersion 
relations and CHPT at next-to-leading order in the chiral expansion. This method 
45| , 0, which does provide a control on the approximations made, would then allow 



for an experimental determination of {a — (3)'"'^ at order E, requiring, however, a 
two-loop evaluation of 77 -^ tt+tt^. 

In Ref. [^, the bound \aTro\ < 35 has been obtained from a study of the e"'"e~ -^ 



7r°7r'^7 reaction. 

Finally, information on the charged pion polarizabilities may be obtained from 
77 -^ 7r+7r~ data in the following manner |^. Both the one-loop expression for the 
transition amplitude and the leading-order expression for a^^ and (3^^ contain the low- 
energy constant Za as the only free parameter. Extracting it from a fit to the cross 



section then determines a-,, and /3^ at this order. The result ^^l ^a = 2.3 ± 1.7 agrees 
within the error with the value /a = 2.7 used in the present work-the corresponding 
numerical values for the leading-order expressions of (a ± /5)*" therefore also agree. 
The cross section in the threshold region is dominated by the Born term contribution 
and is, therefore, rather insensitive to I a ^- This is the main reason for the large 
uncertainty in this determination of /a- 

8.4 Chiral expansion of a^^ and /3.jj- 

The one-loop result is 



[a + pf - 


= , 


a + /5)^ = 


= , 


a-f3f - 


«^"a , 


247r2M,F2 ''■" ' 


a-Pf - 


- 10 


ASn^M^F^ 



(8.12) 

Here we have identified F with the physical value of the pion decay constant, and 
we have used the charged pion mass to evaluate of (a — /J)*". It is straightforward 
to determine from the amplitudes given in the previous section the neutral pion 
polarizabilities to two loops. The numerical results are displayed in table 3. The 
second column contains the contribution at order E~^, and the third to fifth columns 
display the terms of order E. The total values are given in column 6. (The two-loop 
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Table 3: Neutral pion polarizabilities to two loops in units of 10"*^^ . The contri- 
bution due to cliiral logarithms, listed in the fifth column with bracketed numbers, 
is included in the two-loop result quoted in column four. 





0{E-') 


0{E) 


total 


uncertainty 


1 loop 


K 


2 loops 


chiral logs 


{a+^r 


0.00 


1.00 


0.17 


[0.21] 


-1.15 


±0.30 


{a - pr 


-1.01 


-0.58 


-0.31 


[-0.18] 


~ -1.90 


±0.20 


a^o 


-0.50 


0.21 


-0.07 


[0.01] 


~ -0.35 


±0.10 


/3.0 


0.50 


0.79 


0.24 


[0.20] 


~1.50 


±0.20 



contribution (a±/5)^ = 0.18 reported earlier [^ and quoted in Ref. [^ corresponds 
to slightly different values of /i and l2-) Finally, our estimate of the errors is shown in 
the last column. These are obtained in the same manner as the ones for the couplings 
/tj. and /i^ in Eq. ( [6.151 ). We have not considered correlations in these uncertainties, 
which do also not contain effects from higher orders in the quark mass expansion. 

The contribution from the chiral logarithms present in the low-energy expansion 
of the polarizabilities deserves a comment. As we discussed earlier, the In^ M^/ ij? 
terms occur in a particular combination which is dictated by the general structure 
of the renormalized amplitude. 



(« ± /3)1 






\YiMl/^^[\nMl/^? + 2l2] . (8.13) 



^ (167r2F2)2 

Here the ellipsis denotes further single logarithms and terms of order M^, and C± 
are Clebsch-Gordan coefficients. These terms are potentially very large. 



-1.14- 10"^fm^ 



(8.14) 



Mp. It turns out that C+ is small, whereas C_ even vanishes. We have listed 



at /x 

the sum of the In^ M"^/ fj? and In M"^/ jji'^ terms at the scale /x = Mp in the fifth column 

of table 3-these contributions are included in the two-loop result quoted in column 

four. 



The low-energy constants determined in [^ give for the contributions from /ij_ 

(0.0,0.72) vector-mesons (nonet) 
(0.0,0.50) chiral quark model 



(«7rO,/?7ro) 



8.15) 
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The corresponding entries in column 3 of table 3 are slightly different than the ones 



from vector-exchange in Eq. ( 8.15 ), because we do not use the nonet-assumption 
here and include in addition axial-vector exchange. The slight discrepancy with the 
chiral quark model prediction is not serious, because the systematic uncertainties in 
that framework are very difficult to assess. 

Turning now to a comparison with the data, we note that the two-loop result 
for {a — (3)^ agrees within the error bars with the value found by Kaloshin and 
Serebryakov [l^. As for the charged pion case, the complete expression at order E is 
not yet available. The chiral logarithms which occur at this order in the low-energy 
expansion can in principle contribute substantially also here. Therefore, to compare 
the chiral prediction with the data, a full two-loop calculation is required |^ 



9 Comparison with dispersion relations 

In this section we compare in some detail the chiral expansion with the dispersive 
calculation carried out by Donoghue and Holstein [13|. 

Consider the S-wave amplitude 

F{s) = ^jdnH^^{s,t). (9.1) 

We find from the two-loop representation given above 

2 



nCHPT 



sFt 



G^{s) [2F^^{s - Ml) + (.2 - Mt)J{s) + C(s, /,)} 



+ 24^^^-^^^)^^»-l6y> + ^- + ^- 



°F = 


= (16.2F|)2[^^^^' + ^^^ 


/l = 


^ ^^-^{-^^-i 


/2 = 


^ ^^h\^^-^ ' 


/ = 


, M2 
= In , 

/i2 



(9.2) 

where C{s,li) is displayed in Eq. (|7^ ), and where Ap is the 5- wave contribution 
from A^^. In the region 2Mn- < E < 400 MeV, the polynomial Pp contributes 
very little to the amplitude. For the comparison with the dispersive approach in this 
region we therefore drop it, together with Ap- 
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In the simplest version of their analysis, Donoghue and Holstein write 

T-lDISP '^^Tt[S) 



3sF^ 



{{2s - M'^)D,\s) + {s- 2M'^)D,\s)} 



where Dj^ is the Omnes function 

^' ^ 1-kiS- t?n67r J(s) ' ^^'^^ 

with 

1 , 1 



k 







25M2 ' 30M2 



^° " 32nF^ ' 2 - 32vrF| " ^^"^^ 

Expanding D^^ ^^^ keeping terms of the same order as in i^cHPT^ ^^ ^^^ ^^^ count 
/c/ as order E^) 

2 



-GAs) {2F'As - M^2) + is' - Mt)J{s) + C^-} 



/a f , , ,ox ^, X 2 



+ ^^ {(^ - M^^) J(.) + -Fls{h -k,)] + 0{E') , 

2F^s 

C^°"" = ^[fco(25-M2) + fc2(s-2M^2)]. (9.6) 



The two representations ( p.2| ) and (p.6| ) give very similar cross sections up to E' ~ 400 
MeV. This is at first surprising, because the polynomial C(s, /j) in the chiral repre- 
sentation ( |9.2| ) contains rescattering effects which are algebraically quite different 
from C^^^^ (e.g., the leading terms ^ s^ differ by more than a factor of 3). The 
polynomial multiplying /a is also different in the two representations. The combined 
effect of these two differencies is that the S-wave amplitude ( p.6| ) agrees numerically 
quite well with (^]2|). [Notice that Fig. 3 in Ref. |13[ which displays the cross section 



according to Eq. (|973|) is not correct ||5l 



Donoghue and Holstein then refine their representation ( |9.3D by adding contribu- 
tions from resonance exchange. Their final result for the cross section agrees very 
well with our representation below E = 400 MeV, see Fig. 11. There we display 
with a solid line the two-loop result. The dashed line is the result of Donoghue and 
Holstein (Fig. 2 in [0). The two representations differ in the threshold region, 
because M.„ is identified with the charged pion mass by these authors. 
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There are differences in the two representations, though. First, in the dispersive 
method, higher order terms are partially summed up. We consider the fact that the 
cross sections agree as an indication that yet higher orders in the chiral expansion 
do not affect the amplitude in the threshold region very much. Secondly, CHPT 
reveals that the amplitude contains chiral logarithms, generated by pion loops. All 
of these effects are not incorporated in the dispersive analysis of Ref. ||13|]. To 
illustrate, consider the amplitude F at the Compton threshold, where it determines 
the difference of the electric and magnetic polarizabilities, 

{a-(3r = ^m. (9.7) 

Numerically, the chiral logarithms amount to a 18% correction to the leading-order 
term (a — /5)^ = —1.01, see table 3. The result {a — [3)^ = —1.76 quoted in Ref. |13 



corresponds to the one-loop contribution and to vector exchange alone and therefore 
differs from our value {a — (3)^ ^ —1.90. In fi^o, axial-vectors do not contribute. 
The logarithms amount to 0.20 in the final result (3t^o ^ 1.50 which differs by 20% 



from the value /^jrO = 1.26 in [13 
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10 Summary and conclusion 

1. At leading order in the chiral expansion, the amphtude for 77 -^ tt^tt^ is 
generated by one-loop graphs |l|, If. In the case of SU{2)l x SU{2)r x t/(l) 
considered here, it involves the pion decay constant and the pion mass as the 
only parameters. The corresponding cross section deviates from the data and 
from dispersive calculations already near threshold. 

2. The neglected terms in this calculation are related to tttt final-state interactions 
and to three new low-energy constants /i!j_ and /i^ which occur at order E^ in 
the effective action. 

3. To investigate these corrections, we have evaluated the next-to-leading or- 
der terms in the chiral expansion (two-loop diagrams) and have estimated 
the new couplings in the standard manner 0, ^ by resonance saturation 
(J^c^0++,1— ,l+-,2++). 

4. The improved cross section agrees rather well with the data and with dispersion 
theoretic calculations at and also substantially above the threshold region, see 
Fig. 5 and Fig. 11. The enhancement in the cross section is mainly due to tttt 
rescattering and to the renormalization of the pion decay constant. 

5. The two-loop corrections are not unduly large-their size is similar to the cor- 
responding next-to-leading order correction in the isospin zero vrvr scattering 



amplitude [ffl and in the scalar form factor of the pion [21 



6. The couplings /i!j_ and /i^ contribute with a negligible amount below E = 400 
MeV [Q-[^. Above this energy, the inherent uncertainty in /i^ becomes more 
important (Fig. 9). The influence of /i!j_ is quite small also in the region 400 
MeV< E < 600 MeV. 

7. The amplitude for the crossed reaction 77?° -^ 77?° is small at the threshold 
E^T, = M^. As a result of this, the one- loop representation is substantially 
distorted by the next-to-leading order terms, although there are no final-state 
interactions in this case. The dominant effect is due to h[_ (omega exchange 
'jTT^ ^ uj ^ 77r° in the language of resonance saturation 



8. The quark mass expansion of the pion polarizabilities a-j^o and jSt^o contains 
chiral logarithms ~ M^ In M^ and ~ Mt^ In Mt^ which contribute substantially 
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to 0.^,0 ±/5^o, although their effect is suppressed by small Clebsch-Gordan coef- 
ficients. The effect of the low-energy constants h^^ on the value of 0.^,0 =F /^ttO is 
large. It will presumably be difficult to extract these couplings from low-energy 
77 -^ 7r°7r° data alone and to determine in this manner the polarizabilities at 



two-loop order [|13 



9. The DAFNE facility |T^, [T^ will have the opportunity to test the chiral predic- 
tions at next-to-leading order in much more detail than is possible with present 
data. 
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A Decomposition of V^^ 

Here we briefly discuss the correlator 

V^, = tJ dxe-^^'^'^+'^y'^ < 7r°(pi)7r°(p2)out | Tj,{x)Uy) | > . (A.l) 

We decompose V^"^ into Lorentz and parity invariant amplitudes 

V^j,u = Vog^^ + ViQi^qi^ + V2qi^,q2u + V^Qi^A^ 

+V4qiuq2fM + V^qi^Af, + V6q2^,q2u + V7q2^,A^ 

A^ = {Pl - P2)t, 

Vi = V,{s,v) ; z = 0,...,9 , 

s = (gi + ^2)^ , t = {pi- qif , M = (p2 - qif , 
V = t-u. (A.2) 

From Bose symmetry 

V;,(A, gi, ga) = V;,(-A, q,, ga) = K^(A, ga, gi) (A.3) 

we find 

V,{s,u) = Vi{s,-iy) ■ 2 = 0,1,2,4,6,9 

V,{s,i^) = -V,{s,-u) ; 1 = 3,5,7,8 (A.4) 

and 

Ve = V,,V5 = -Vr , Vs = -V, . (A.5) 



The Ward identity 

gives with qf = 



q'^V^. = (A.6) 

2^0 + sV^ -uV5 = , 
sVi + iyV3 = , 
sVj -iyVg = 0. (A.7) 

The second Ward identity q2Vfj,u = is then automatically satisfied by Bose symme- 
try. We are left with four independent form factors which we take to be 

A = -V4 , B = Vg/2s , C = V2 , D = Vs/s . (A.8) 

Insertion into (|A.2| ) gives the decomposition Eq. ( p^) in the text. 
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B 77 ^ TTTT to one loop 

For convenience, we collect here the one-loop expressions for the amplitudes 77 



7r°7r°, TT'^TT 



Bi- 



B.l 77 -^ TT+TT 

The matrix element for 

7(^1)7(92) -^ 7r+(pi)7r"(p2) (B.l) 

is given by 

< 7r+(pi)7r-(p2)out I 7(gi)7(g2)in >= t{2ny6\Pf - P,)T'' , (B.2) 

with 

V^^ = i /" rfxe-*(''^^+«^^) < TT+TT-out I Tj^{x)3^{y) \ > 

= A^Ti^, + E^Ts^, + C^'Ts^, + D^T4^, . (B.3) 

The tensors Tj^^ are defined in ( p.4|) . The hehcity amplitudes are 

s 

The low-energy expansion of the amplitudes A'~^, B'-'' reads [^ in SU{2)l x SU{2)ji x 
f/(l) with the Condon-Shortley phase convention 



The loop-function G{s) is discussed in appendix 0. We do not split the result into 
A2,A4 etc., because the propagators contain the physical pion mass at one-loop 
order-this would make the splitting rather useless. 
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B.2 77 -^ ttOtt^ 

The leading term is generated by one-loop diagrams alone-there is no contribution 
from £4. The result is [|], § 

54 = . (B.6) 

CHPT thus predicts the cross section at this order in the energy expansion in terms 
of the two parameters F and M^. The amplitude is purely S-wave. 

In order to compare the prediction (|B.6|) with the data, we identify F{M) with 
the physical pion decay constant F^ (physical pion mass M^), as this induces only 
changes of higher order. The result is shown in Fig. 5, where we display the cross 
section a{s] \ cos 6 \< 0.8) according to Eq. ( p.6|) with a dashed line, together with 
the Crystal Ball data P] as a function of the center-of-mass energy E = ^/s. The 
cross section is below the data for E < 400 MeV where the low-energy expansion 
can be trusted most. It also differs by a similar amount from dispersion theoretic 
calculations |[lO|-ll^. An example (Fig. 23 in |12|) is shown as dash-dotted lines in 



the figure. The solid line is the two-loop result. 

The amplitude ( p.6|) has the peculiar property that its dispersive representation 
needs a subtraction, although the absorptive part vanishes at high energy sufficiently 
fast to generate a convergent unsubtracted dispersion integral. 



lmA = 0[ ,s^oo. (B.7) 




Finally, we note that the leading term ( [B.6| ) approaches a constant in the chiral 
limit, 

^(^) = -i;^ + 0(^') '^^0 '^^0- (B-8) 



C Loop-integrals 

1. The loop-integral G{s) is 
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G is analytic in the complex s - plane, cut along the positive real axis for Re s > 4M^. 
At small s, 



n! 



G{s) ig7r2^^(^M2J (n + l)(2n+l)! ■ 



The absorptive part is 

ImG'(s) 



M^ r 1 + (7 
m 



STT I 1 — C 



s>4M2 , 



a = Jl-AM^/s. 



Use of 



U2{y) + Li2 



-y 



\n\l-y) 



U,{y) = -r^ln(l-x) 
Jo X 



ry dx 



gives 



16712^(5) 



1 +^ ln{^+^^ 



4M^ 



arctg' 



2( s 



^iM'^-s' ' 



1 +Mlln^£^ 



4M2 < s 
< s < 4M2 
s < 0. 



In the text we also need 

G{s) = G{s)-sG'{Q). 

2. The loop-integral J{s) is 

1 /■! s 



IQtt^ 



(C.2) 



(C.3) 



(C.4) 



(C.5) 



(C.6) 



(C.7) 



J is analytic in the complex s - plane, cut along the positive real axis for Re s > 4M^. 
At small s, 



J(s) 



-t(-) 



n! 



167r2;^i VmV n(2n+l)! ' 



The absorptive part is 
Explicitly, 

IGtt^ J(s) = 



ImJ(s) 



a 



IGtt 



s > AM^ . 



(C.8) 
(C.9) 



2-2(M^;- 
^lnS + 2 



2 arctg ( 



4M2- 



4M2 <s 

0<s<4M2 (CIO) 

s < 0. 
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In the text we also need 

J{s) = J{s)-sJ'{0). (C.ll) 

3. The loop-function H is defined in terms of G and J, 

H{s) = {s- 10M^)J{s) + 6M^G{s) . (C.12) 

D Low-energy constants from resonance satura- 
tion 

Here we give details of the calculation needed to estimate the renormalized cou- 
plings a'[,a2 and 6*". We consider the exchange of scalar, (axial-) vector and tensor 



mesons with mass Mji < 1.2 GeV and follow the procedure outlined in Ref. ||30 |. 
The contributions of the vector and tensor mesons are evaluated in the framework 
of SU{2)l X SU{2)ii X U{1). In order to overcome the limitations about the experi- 
mental information presently available on the 1^~ and O"*"^ multiplets we will work 
in SU{3)l X SU{3)r at large Nc- 

D.l Vector and tensor mesons {J^^ = 1 ,2++) 
D.1.1 The lagrangian 



We set 



y.(^-)^\ ? 2\- \7' , (D.l) 



^ i^'-i . 1"='^. 



and have for the kinetic part 



CUV, T)= -\ T.iV.uV^" - 2M^V,V^) 
^ V 



2 M!^ -tpo- 



(D.2) 
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where 



V^, = D^V,-D,V^ , 



(D.3) 

T^ = -\u^ [d^-ir^]u + u[d^-il^]v)^ . 

Furthermore, T^'^ = T^^ denotes the spin-2 field for /2(1270) with J^'-' = 2+'*', and 



:d.4) 



-gP^d^'d" + g'^'d^d" - -{g'^^d^'dP + g'^Pd^'d^ + g^^d^^dP + g^'P^''^''] 



The propagator for T^'^ is obtained in the standard manner by exposing the 
system to an external perturbation, 



We find 



with 



£ = £,„(T)+/%, . (D.5) 

p = -iP P +P P ) P P (D.6) 

p" ^p"^ Ml ' 

DP^-^p''G,^r\^) = lig'^'^g"^ + 5'^V^)<5^(a:) . (D.7) 



Now consider the couplings of V, T to pions and to photons, linear in the reso- 
nance fields. Since we are interested in terms of order E^ in the effective action, it 
suffices to construct interactions which are at most of order E^. We set /+^ = 2eQF^^ 
where F^^ is the photon field, and take 

CUV^T) = ee^,p„FP''Y.{C^y{VP{u'^,Q}) + CUVPu'^){Q)} 

(D.8) 
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where B^'' {Qit'^) is the energy-momentum tensor of the pion (photon) field, 



er = y(«'«') - ^g'^ii^"^^) + (x+)} 



Qf,u ^ pt^^p^^ + -g^"'FP''F. 



(D.9) 



pa 



The couphng V^ — > 7r7 ( T ^ vrvr, 77) has been considered in ref. [^ ([0) for 
the case of nonet fields and (Q) = 0. The interaction (p.8| ) for the spin-2 field differs 
from the one proposed in [^. In particular, our amplitude for 77 -^ tt^tt^is smooth 
at large momenta and purely D-wave also off the /2-resonance, see below. 



D.1.2 The amplitudes 

We find for the amplitudes from vector exchange 
Av = Cv 



s-4(t + M2) s-4(m + M2) 



B. 



C. 



M^-t 



1 



M2 



u 



(D.IO) 



M^-t Ml 



u 



; V = p,uj, 



where 0, ^ ^ 0, ^ 



Cv = -Mi 



T{V -^ TT^) 



(D.ll) 



From the published [^ widths we obtain 

C^ = 0.67 GeV"^ , Cp = 0.12 GeV"^ , C^ = 0.2 ■ lO^^ GeV"^ . (D.12) 

The calculation of the p, u correction to the 77 -^ -k^-k^ scattering amplitude has 



been extended recently to include both photons off shell [Q . 
Tensor exchange gives0 

At + 2{AMl - s)Bt = , 



Br, 



1 r^l f^-K 

AM^-s 



(D.13) 



^ The interaction ( |D.§| ) generates tadpole diagrams where T^ disappears in the vacuum. These 
graphs give rise to additional contributions to the amplitudes, which however are of higher order in 
the energy expansion. We omit them here. 
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with 



r(T ^ 77) 



e^ClYM^ 



yrp 



SOvr 



O^On 



r(r ^ 71%' 



{C^fM^ / 4M2 



, 5 
2\ 2 



9607r 



Ml 



From the measured widths r(/2 -^ 77), r(/2 ^ tt^tt*^) = |r(/2 

\Cl\ = 0.19 GeV"^ , 



|CJ| = 9.2 GeV-i . 



(D.14) 



tttt) ||53|, we find 



(D.15) 



D.2 Axial-vector and scalar mesons (J^^' = 1+ ,0^^) 
D.2.1 The lagrangian 

In this paragraph, A denotes a 3 x 3 matrix. In particular, 



5. 



u. 



V2 



F ^- ' 








V2^ VE 


7r+ 


K+ 


\ 


7r^ 


V2^ Vq 


K' 




K- 


K^ 


u 


1 



(D.16) 



and Q = |diag(2, —1, —1). (We do not include 77-77' mixing, as this is of higher order 
in the quark mass expansion.) Furthermore, B^ stands for the axial-vector nonet. 



^.(1^~) 



V2 



BIX 



-^^M-1 
v^ ^ 



(D.17) 



(D.18) 



and similarly for the scalar nonet S'(0^"'"). The kinetic terms are 

CUS,B) = ]^{D^SD^S-MlS^) 

-^{B^.B^^" -2MlB^B^) , 

where Mb denotes the common nonet mass, and the covariant derivatives are the 
SU{?>) version of ( p.3| ). The couplings to pions and to photons are |^, ^, ^ 

fD 19) 
+ C^{Sx+)+eCBF,,{B>^{Q,u''}) . ^ ' ^ 
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D.2.2 The amplitudes 

We find[^ for the contribution from the scalar exchange [2^ 

As 






9F^{MI - s) 



Bs = 



with 



r(ao -> 77) 



72rr 



The couplings C^ and C^ have been determined in |30 



\Ci\ =3.2- 10"' GeV 



(D.20) 



(D.21) 



I C™ I =4.2-10-" GeV , 



(D.22) 



C^C^ > . 

The quantity \Cg\ is difficult to estimate. Here we content ourselves with a rough 
estimate by relating it to the decay Oq — > 77 via 

r 



r(ao -> 77) = r- 



(D.23) 



r(ao -^ -n^ri) 

where r = 0.24 keV [Q with a sizeable error. We furthermore assume that the decay 
Qq — > TT^r] accounts for all of r^^t 



Hence we obtain 



r(ao -^ 1x^7]) 



\Cl\ = 8.2 ■ 10-2 GeV"^ . 



(D.24) 
(D.25) 



For the contribution from the axial nonet exchange |50, 38 1 we find 



Ab = Cb 



s + A{t- Ml) s + A{u- Ml) 



B 



C 



B 



B 



Ml-t 



1 1 

+ 



Ml-u 



M^-t Ml-u 



(D.26) 



^° The interaction ( D.igj ) generates tadpole diagrams where {S) disappears in the vacuum. These 
graphs give rise to additional contributions to the amplitudes, which however are of higher order in 
the energy expansion. We omit them here. 
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with 

_30 3 r(6i^7r+7) .^271 



or 

-2 



Cs = 0.53 GeV"' . (D.28) 

D.3 Expressions at low energies 

At low energies, the above contributions from V, B, S and T sum up to 
Aq = aiMl + 0,28 , 






(D.29) 



Cy _ r, ^g I / 20 |CgCg| 1 |CyCy| 
^^^^^^^ M2 Ml \9 F^MI 2 M| 



«2-6 J2 172" ~ 2— ^± — /2 + 






v^.,<^M'v Ml A M^ ■ 

In table 2 we list the contributions from the individual resonances to the dimension- 
less parameters (167r^F^)^(ai; a2; b). We used the following values for the resonance 
masses: M^ = 782 MeV, Mp = 768 MeV, M^ = 1020 MeV, Mb = 1232 MeV, 
Ms = 983 MeV, Mr = 1275 MeV. 
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Figure captions 



Fig. 1. Mandelstam plane showing the three related physical regions, s-channel: 

77 -—>■ TTTT, t-and u-channel: 77? —>■ 'jir. We indicate the threshold for 77 -^ 7ni{'~f7T -^ 
77r) by A (P). The shaded lines at s,t,u = 4M^ indicate the presence of a branch- 
point in the amplitude, generated by two-pion intermediate states. 

Fig. 2. The graphs at order E^ in the chiral expansion. These graphs correspond 
to the terms in Eq. ( p.22|) . The framed symbols I stand for the vertices in Cj. We 



indicate with the solid-dashed lines the propagator D2 ^■ 

Fig. 3. One class of Feynman diagrams which contribute to 77 -^ tt^tt^. The 
dashed (solid) lines stand for neutral (charged) pions. The four-point function on 
the right-hand side is the elastic vrvr scattering amplitude at one-loop accuracy in 
d-dimensions (with two legs off-shell), and the symbol d^l stands for integration over 
internal lines with weight (|5.1|) . 

Fig. 4. Graphs which contribute to Aa,b in the unitary part Ua.b in Eqs. 
73|,[7.8|). In Fig. 4a, the dash-dotted lines surround the diagram which we rep- 



resent in the dispersive manner (|5.3| ). The graph Fig. 4b is called "acnode graph" 

m- 

Fig. 5. The 77 -^ 7r°7r° cross section a{\ cos6\ < Z) as a function of the center-of- 
mass energy E at Z = 0.8, together with the data from the Crystal Ball experiment 
IP]. The solid line is the full two-loop result, and the dashed line results from the 
one-loop calculation |jl], 0. The band denoted by the dash-dotted lines is the result 
of the dispersive calculation by Pennington (Fig. 23 in |ll2|). 

Fig. 6. The dependence of the 77 — » tt^tt^ cross section cr(| cos6\ < Z) on the 
constants k, at Z = 0.8. The sohd hne denotes the two-loop result with the standard 
values for /j displayed in table 1 (without contributions from resonance exchange, 
and with Aa,b = 0), whereas the dashed line is evaluated at k = 0. The dash-dotted 
line has /i^s = and the other /j at their standard values. 

Fig. 7. Real and imaginary part of the amplitudes ±10'^M^H^± at t = u. The 
solid line is for 10^M^if++. It incorporates all contributions except Aa,b- The 
dashed line is the same amplitude for the one-loop case, and the dash-dotted line is 
for — 10^M|if+_ with the same input as the solid line. Finally, the crosses refer to 
the center-of-mass energy of the tt^tt^ system in 100 MeV steps. 
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Fig. 8. The 77 —>■ tt^tt^ amplitude as a function of s/M^ ai t = u. For s < AM^ 
the quantity shown is 10"^ M^H++ and for s > AMI we display lO^M^Iif+^l. The 
solid line is the result of the two-loop calculation and the dashed line is the one-loop 
result. The symbols (0,+) refer to the work of Pennington |1^: o{+) is from Fig. 
19 (23) in that article. 

Fig. 9. The uncertainty in the value of the 77 -^ tt^tt^ cross section o"(| cos 0| < Z) 
at Z = 0.8, evaluated from the two-loop amplitude at ^a,b = 0. The data are from 
the Crystal Ball experiment [^j. The dashed lines embrace the region generated by 
assigning all possible combinations of signs to the systematic errors in the couplings 
/ij_ and /i^ according to Eq. ( |6.15| ). The dash-dotted line corresponds to the central 
values in (|6.15|) . Above E ~ 400 MeV, the major part of the uncertainty in the cross 
section is generated by the error in /i^. 

Fig. 10. The Compton cross section 77r° — > 771° as a function of the center- 
of-mass energy E^^^. The solid line is the result of the two-loop calculation and 
the dashed line is the one-loop result. The dash-dotted line refers to the two- loop 
calculation at /i!J_ = /i^ = 0, and the dotted line has iJ+_ = in a full two- loop 
calculation. 

Fig. 11. The 77 -^ -k^-k^ cross section o"(| cos^| < Z) as a function of the center- 
of-mass energy at Z = 0.8, with the data from the Crystal Ball experiment. The 
solid line is the two- loop result, whereas the dashed line is taken from the dispersive 
analysis of Donoghue and Holstein (Fig. 2 in |T^). 
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